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ABSTRACT 

We present new near infrared spectroscopic measurements of the H/3 region for a sample of 29 luminous high 
redshift quasars. We have measured the width of H/3 in those sources, and added archival H/3 width mea- 
surements, to create a sample of 92 active galactic nuclei (AGNs) for which H/3 width and rest-frame UV 
measurements of N V A 1240 and C IV A 1549 emission-lines are available. Our sample spans six orders of 
magnitude in luminosity and includes 31 radio-loud AGNs. It also includes 10 narrow-line Seyfert 1 galaxies 
and one broad absorption-line quasar. We find that metallicity, indicated by the N V/C IV line ratio, is primarily 
correlated with accretion rate, which is a function of luminosity and H/3 line-width. This may imply an intimate 
relation between starburst, responsible for the metal enrichment of the nuclear gas, and AGN fueling, repre- 
sented by the accretion rate. The correlation of metallicity with luminosity, or black hole (BH) mass, is weaker 
in contrast with recent results which were based on measurements of the width of C IV. We argue that using 
C IV as a proxy to H/3 in estimating Mbh might be problematic and lead to spurious BH mass and accretion 
rate estimates in individual sources. We discuss the potential implications of our new result in the framework 
of the starburst-AGN connection and theories of BH growth. 



Subject headings: galaxies: active - galaxies: nuclei 
abundances - galaxies: starburst 

L INTRODUCTION 

The determination of chemical abundances in active galac- 
tic nuclei (AGNs) provides a powerful tool to probe star form- 
ing activity in galactic nuclei at present, and in the early 
Universe. Early metal abundance studies in AGNs utilized 
weak and broad intercombination lines, such as N III] A 1750, 
N IV] A1486, and O III] A1663, to determine metal abun- 
dances in the broad-line region (BLR) gas (Shields 1976). 
These early studies have found that BLR metallicity is about 
solar, with some luminous quasars reaching super-solar val- 
ues. Hamann & Ferland (1993; hereafter HF93) found that 
emission-line ratios, such as N V A1240/C IV A 1549 (here- 
after N V/C IV) can trace the BLR metallicity across a wide 
range of luminosity and BLR density (see also Hamann et 
al. 2002 for a more detailed model). HF93 also found that 
BLR metallicity, determined in this way, correlates with lu- 
minosity, in what they proposed as the metallicity-luminosity 
(Z-L) relationship in AGNs. Given the well known rela- 
tionship between AGN luminosity and black hole (BH) mass 
(Mbh), HF93 suggested that the Z-L relationship can natu- 
rally lead to a Z-Mbh dependence, in analogy with the mass- 
metallicity relation observed in some elliptical galaxies (e.g., 
Trager et al. 2000). 

Shemmer & Netzer (2002; hereafter SN02) have shown that 
once narrow-line Seyfert 1 galaxies (NLS 1 s) are introduced to 
the Z-L diagram, they deviate significantly from the Z-Lxt- 
lation by exhibiting high N V/C IV at low luminosity. A possi- 
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ble explanation is that N V/C IV is not an adequate metallicity 
indicator for NLS Is (and perhaps other AGNs). Otherwise, 
this implies that at least at low-luminosity, the Z-L relation is 
more complex and cannot be a simple two-parameter depen- 
dence. According to SN02, BLR metallicity also depends on 
the width of H/3, which is perhaps the best accretion rate (in 
terms of the Eddington ratio: LBoi/i'Edd) indicator, and that 
this dependence may prevail also at the high-luminosity end. 
To test this hypothesis, we observed the near infrared (IR) 
spectrum of a sample of 29 luminous 2 < z < 3.5 quasars, 
with available N V/C IV data, in order to measure the width 
of their H/3 lines, and thus deduce the accretion rate in these 
sources. 

This paper is part of a two-paper project describing near- 
IR spectroscopic measurements of luminous high-redshift 
quasars aimed at estimating iBoi/i'Edd and measuring the 
emission line properties in such objects. In this (first) paper 
we describe the quasars' selection, their observations and data 
reduction (§|2ji. A description of the data analysis is given in 
§ E In § |4] we summarize our main results concerning the 
broad emission lines, and in §|5lwe discuss their potential im- 
plications. A summary of our main conclusions is given in 
§ |6l A following paper (Netzer et al. 2004; hereafter Paper 
II) addresses the complementary issue of the [O III] and Fe II 
lines in those sources. 

2. SAMPLE SELECTION, OBSERVATIONS, AND DATA 
REDUCTION 

We have selected a sample of 29 luminous quasars for our 
study upon the following criteria: 

1. Luminous (L ^ 10^^ erg s~') sources with expected H 
magnitudes ^ 17 to allow high-S/N spectra on medium 
size telescopes. 

2. Each source has an archived or published rest-frame 
UV spectrum that includes the N V and C IV emission 
lines, which are not subject to severe absorption. 
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Table 1. Observations Log 



Quasar Name 


RA (2000) 


DEC (2000) 


z" 


Z.y.' 


H^ 


Observatory 


Date 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


2QZJ00 122 1.1 -283630 


00 12 21.1 


-28 36 30 


2.327 


2.339 


16.6 


AAT 


2003 Jul 17 


2QZJ002830.4-281706 


00 28 30.4 


-28 17 06 


2.406 


2.401 


15.9 


AAT 


2003 Jul 16 


UM667 


00 47 50.1 


-03 25 31 


3.122 


3.132 




TNG 


2003 Aug 18 


LBQS 0109+0213 


01 12 16.9 


+02 29 48 


2.343 


2.349 


15.3 


TNG 


2002 Nov 8 


[HB89] 0123+257'' 


0126 42.8 


+25 59 01 


2.358 


2.369 


15.9 


TNG 


2002 Aug 27 


2QZ J023805. 8-274337 


02 38 05.8 


-27 43 37 


2.452 


2.471 


16.0 


AAT 


2003 Jul 15 


SDSS J024933.42-083454.4 


02 49 33.4 


-08 34 54 


2.491 


2.491 


16.5 


AAT 


2003 Jul 16 


[HB89] 0329-385'' 


03 3106.3 


-38 24 05 


2.423 


2.435 


15.6 


AAT 


2003 Jul 14 


[HB89] 0504+030'' 


05 07 36.4 


+03 07 52 


2.463 


2.473 




TNG 


2002 Nov 10 


SDSSJ100428.43+001825.6 


10 04 28.4 


+00 18 26 


3.046 


3.046 




TNG 


2003 Nov 13 


TON 618'' 


12 28 24.9 


+31 28 38 


2.219 


2.226 


13.9 


TNG 


2003 Jun 4 


[HB89] 1246-057 


12 49 13.9 


-05 59 19 


2.236 


2.240 


14.3 


TNG 


2003 Jun 6 


[HB89] 1318-113 


13 21 09.4 


-1139 32 


2.308 


2.306 


15.1 


TNG 


2003 Jun 5 


[HB89] 1346-036 


13 48 44.1 


-03 53 25 


2.344 


2.370 


15.1 


TNG 


2002 Apr 5 


SDSS J135445.66+002050.2 


13 54 45.7 


+00 20 50 


2.511 


2.531 




TNG 


2004 Apr 1 


UM629 


14 03 23.4 


-00 06 07 


2.462 


2.460 


16.0 


AAT 


2003 Jul 16 


UM 632'' 


14 04 45.9 


-01 30 22 


2.518 


2.517 


16.1 


TNG 


2003 Jun 5 


UM642 


14 10 26.4 


-00 50 09 


2.372 


2.361 


16.6 


TNG 


2003 Jun 4 


UM 645'' 


14 1123.5 


+00 42 53 


2.269 


2.257 




TNG 


2003 Jun 7 


SBS 1425+606 


14 26 56.1 


+60 25 50 


3.160 


3.202 


14.5 


TNG 


2003 Aug 30 


SDSS J170102.18+612301.0 


17 0102.2 


+61 23 01 


2.293 


2.301 


16.4 


TNG 


2003 Jun 4 


SDSS J173352.22+540030.5 


17 33 52.2 


+54 00 30 


3.425 


3.428 


15.7 


TNG 


2003 Aug 30 


[HB89] 2126-158 


21 29 12.2 


-15 38 41 


3.268 


3.282 


14.9 


TNG 


2003 Oct 15 


[HB89] 2132+014'' 


21 35 10.6 


+01 39 31 


3.194 


3.199 




TNG 


2003 Sep 5 


2QZJ221814.4-300306 


22 18 14.4 


-30 03 06 


2.384 


2.389 


16.0 


AAT 


2003 Jul 16 


2QZ J222006.7-280324 


22 20 06.7 


-28 03 24 


2.406 


2.414 


14.3 


AAT 


2003 Jul 15 


[HB89] 2254+024'' 


22 57 17.5 


+02 43 18 


2.081 


2.083 


15.9 


TNG 


2003 Jun 16 


2QZJ231456.8-280102 


23 14 56.8 


-28 01 02 


2.392 


2.400 


16.6 


AAT 


2003 Jul 17 


2QZJ234510.3-293155 


23 45 10.3 


-29 3155 


2.360 


2.382 


16.5 


AAT 


2003 Jul 17 



"Redshift taken from the NASA/IPAC Extragalactic Database (NED; http://nedwww.ipac.caltech.edu/ i, except for boldface values that were taken from the 
2QZ archive. 

''Systemic redshift determined from the peak of the [O III] A5007 emission line. In cases where no [O III] emission was detected, H/3 was used instead (bold 
face). 

'^Magnitudes in bold face are 2MASS measurements, while others were obtained by photometric measurements using field comparison stars, whose magnitudes 
were taken from the 2MASS archive. Six sources do not have 2MASS magnitudes (or H images with comparison stars having 2MASS magnitudes), and their 
rest-frame optical fluxes were obtained from extrapolation of the rest-frame A 1450 flux. 

''Radio loud quasar. 



3. H/3 is located in the observable IR bands, where it is 
least affected by atmospheric absorption. 

According to these criteria all of our quasars lie in a 2 < 
z < 3.5 redshift range, with H/3 located in either the H or K 
bands. Observations were done at the Anglo-Australian Tele- 
scope (AAT) in Australia and at Telescopio Nazionale Galileo 
(TNG) in Spain. An observations log and quasar details ap- 
pear in Tabled 

The observations at TNG were done with the Near Infrared 
Camera Spectrometer (NICS). Eight quasars were observed 
with a l"-wide slit, and eleven with a 0.5"-wide slit. During 
the observations the telescope was nodded along the slit and 
each spectrum was taken at 2-6 different positions along the 
spatial axis in order to perform the primary background sub- 
traction. We used the HK grism to obtain the 1.4-2.5 /im 
wavelength range, which provides an almost uniform disper- 
sion of 11.2 A/pix throughout this range. This results in a 
resolving power of R - 500 (1" sUt) or R - 1000 (0.5" sHt), 
slightly increasing towards the red. Each night, spectra of one 
or two standard O stars were taken as close as possible in time 
and in airmass to the quasars to allow the removal of telluric 
absorption features from the quasars' spectra. 

The spectroscopic observations at AAT were done with the 
second generation Infrared Imager & Spectrograph (IRIS2) 



and a l"-wide slit. During the observations the telescope 
was nodded along the sht in an ABBA sequence, thus each 
spectrum was broken into a set of 4 frames to allow primary 
background subtraction. We used the Sapphire-H Grism to- 
gether with an H filter to obtain the 1.54- 1.90 /im wave- 
length range^, with a dispersion of 3.4 A/pix, and a resolution 
of R ^ 2500. Standard stars of spectral types O, B, and G 
were observed each night to remove telluric features from the 
quasars' spectra. We also used the imaging mode of IRIS2 
to take H band images of our quasars. H magnitudes for the 
quasars were obtained by relative photometry and comparison 
with H magnitudes of nearby stars taken from the two micron 
all sky survey (2MASS) archive *'. Seven of the quasars also 
had H magnitudes in the 2MASS archive. Non-systematic 
differences of up to ^ 0.5 mag were found between these 
magnitudes and the magnitudes obtained by the relative pho- 
tometry. This may be mainly due to increased photometric 
errors near the edge of the 2MASS detection limit. 

^ Normally, this setup allows a coverage of the 1.485 - 1.781 ^m range, 
i.e., the entire H band, however, a misplacement of the grism in the filter 
wheel during our run, resulted in a short wavelength cutoff. 

* http://www.ipac.caltech.edu/2mass/ 
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Fig. 1 . — IR spectra of 2 < z < 2.5 quasars obtained at AAT. The spectrum in eacli panel is represented by a thin solid line. The model (thick line) is composed 
of a continuum {dotted line), Fe II emission, H/3, and [O III] (long dashed lines). 



The two-dimensional spectra (quasars and standard stars) 
from both observatories were reduced using standard IRAF^ 
routines. All spectra were pair-subtracted to remove most of 
the background, which is swamped by strong telluric features, 
such as OH emission, leaving some residuals perpendicular to 
the spatial axis. For each observatory we used a lamp spec- 
trum as a flat-field image to correct each spectrum for pixel- 
to-pixel variations. Each spectrum was fitted by a low order 
polynomial, and a secondary background subtraction, along 
the dispersion, was carried out. An extraction window of 1" 
(0.5") was used to extract each of the AAT (TNG) spectra, 
respectively. Wavelength calibration was achieved using Xe 
and Ar arc lamps. Finally, all spectra were rebinned to a com- 
mon resolution of ^ 600 km s~\ matching the resolution of 
the TNG spectra taken with the 1" slit. 

The spectral shape of the quasars was recovered by divid- 
ing each spectrum by a spectrum of a standard star, after 

' IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by AURA, Inc., 
under cooperative agreement with the National Science Foundation. 



the stellar absorption features were carefully removed. The 
quasar/star ratio was then multiplied by a stellar model in the 
form of a black body curve, with a temperature corresponding 
to the star's tabulated spectral type. Some of our standard stars 
have spectra that appear in the library of near-IR stellar spec- 
tra of Lan^on & Rocca-Volmerange (1992). This enabled the 
removal of major intrinsic absorption features in our standard 
stars, and also allowed to compare our stellar model to the 
archived data. Flux calibration of the quasars was obtained 
by comparing synthetic H and K magnitudes of the quasars 
(by integrating counts across the dispersion) to the observed 
H magnitudes (see Table[0, and using the flux density of Vega 
as the zero-point flux. 

Our flux calibration may be subject to relatively large un- 
certainties emerging from different atmospheric absorption 
between a quasar and its corresponding standard star. This 
uncertainty increases for the AAT spectra, as those are cutoff 
short of the H band blue end. For those spectra we integrated 
the predicted counts after extrapolating the continuum (fitted 
by a linear function) beyond ^ 1 .55 /im. We also used the 
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Fig. 2. — IR spectra of 2 < z < 3.5 quasars obtained at TNG. Symbols are as in Fig.fTI 
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bj magnitudes of the 2QZ quasars to check for consistency 
with our derived rest-frame optical fluxes, by converting those 
magnitudes to observed 4400A fluxes and extrapolating them 
to 5100A assuming a fiducial quasar continuum function of 
the form: f,y °^ i^"*^^. We found that our H and by magnitudes 
(reflected to the observed H band) were consistent to within 
±0.25 magnitudes. We note that those two magnitudes were 
measured at different epochs, separated by ^ 20 years in time, 
thus constraining the variability amplitude of each quasar to 
~ 25% over a time scale of ^ 6 yr in the rest-frame (although 
there is a scope for variability within this long period of time). 
Luminosities were obtained by using the following cosmolog- 
ical parameters: Q\ = 0.7, Qyn = 0.3, and Hq = 70 km s~' 
Mpc"^ 

In addition to our 29 sources, we added to our sample two 
quasars, LBQS 0256-0019 and LBQS 0302-0000, for which 
rest-frame UV spectra were available from the Sloan Digi- 
tal Sky Survey (SDSS; Stoughton et al. 2002) archive, and 
a high-quality rest-frame optical spectrum including H/3 was 
kindly provided by M. Dietrich (private communication; see 



Dietrich et al. 2002). As a result, our high redshift sample 
consists of 31 luminous quasars. 

Our new sample was complemented by 61 z < 2 AGNs, 
from the original SN02 sample, all with published Hfi data. 
Altogether this sample comprises 92 sources, including 10 
NLSls, with N V, C IV, and H/3 information. One of our 
sources is a broad absorption line quasar ([HB89] 1246-057), 
based on its UV spectrum, and 3 1 of our sources (8 at 2 < z < 
3.5, and the rest at low-z including one NLSl) were identified 
as being radio loud quasars (RLQ), following the discrimina- 
tion between RLQ and radio quiet quasars (Kellerman et al. 
1989), and using optical and radio fluxes from Veron-Cetty & 
Veron (2003). 

3. MEASUREMENTS AND DATA ANALYSIS 

3.1. HP width 

The H/3 region in the final spectra was fitted with a model 
consisting of a linear continuum, the Boroson & Green (1992) 
Fe II emission template, and a multi-Gaussian fit to the emis- 
sion lines. We also fitted the Ha region in the TNG spec- 
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Fig. 3. — IR spectra of 2 < z < 3.5 quasars obtained at TNG (continued). Symbols are as in Fig.l2l 



tra (where present) in the same way, but these data are left 
out of the current analysis. In the first step we fitted a lin- 
ear continuum and one Gaussian to H/3 to get a rough esti- 
mate of its full width at half maximum (FWHM). The next 
step was to broaden the Fe II emission template by con- 
volving it with a Gaussian whose width was identical to 
FWHM(H/?). We then fitted the spectrum with a linear con- 
tinuum, the broadened Fe II template, shifted in flux by a 
scaling factor, and multi-Gaussian components for the emis- 
sion lines. Two Gaussians with differing widths were as- 
signed to H/3 and two Gaussians with identical widths were 
assigned to the two [O III] AA4959,5007 Hnes. The fluxes 
of the [O III] lines were constrained to the theoretical ratio 
I([OlII] A5007)/I([OIII] A4959)= 2.95, and the widths of the 
two Gaussian components of H/3 were free to vary. 

The final spectra, best fit model, and model components, 
are plotted in Figures ^^3] These fits are the ones used for 
FWHM(H/3) listed in Tablegl The final FWHM of the entire 
H/3 profile was measured, as well as the FWHM of the [O III] 
lines. We estimated the uncertainty on the width of the lines 
by slightly shifting the parameters from their best-fit values 



and marking the most extreme values that would still repre- 
sent a reasonable fit to the H/3 complex. In almost all of the 
cases, those uncertainties were smaller than the instrumen- 
tal resolution (~ 600 km s~') for the TNG and AAT spectra, 
with UM 645, [HB89] 2254+024, 2QZ J002830.4-281706, 
and 2QZ J23 1456.8-280302 being the only exceptions (see 
Figures n|3j. 

We have also obtained direct measures of the fluxes and 
FWHMs of H/3 and [O III], by integrating across the Fe II- 
subtracted features, for the flux measurements, and by taking 
the width at half the peak value for each Fe Il-subtracted fea- 
ture, in the FWHM measurements. Both types of measure- 
ments, the best-fit value, and the 'direct' feature method, ap- 
pear in Table |2] and in Paper II. The difference between both 
methods can be regarded as an estimate of the uncertainty on 
those measurements, although these are not pure statistical er- 
rors. This paper presents the H/3 line width for our sources. 
Other line properties, and in particular those of [O III] and 
Fe II, are discussed in Paper II. 

It is expected that each H/? profile in our 2 < z < 3.5 sample 
is a blend of a narrow component, emitted from the narrow- 
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Table 2. Emission Line and Continuum Measurements 



Quasar Name 


Log 


FWHM(H/3) 


Log 


L/LEdd 


N V/C IV 


Reference 




AL;^(5100) 


Best Fit 


Direct 


Mbh 






for 




[ergs-l] 


[kms-'] 


[kms-'] 


[M@] 






UV data 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


2QZJ001221. 1-283630 


46.26 


1915 


2382 


8.71 


1.64 


0.84 


1 


2QZJ002830.4-281706 


46.58 


4833 


5295 


9.71 


0.35 


0.47 


1 


UM667 


46.28 


3135 


2244 


9.15 


0.62 


0.81 


2 


LBQS 0109+0213 


46.80 


5781 


5241 


10.00 


0.30 


0.82 


3 


[HB89] 0123+257 


46.57 


2406 


2763 


9.09 


1.38 


0.77 


2 


2QZ J023805. 8-274337 


46.57 


3437 


3691 


9.40 


0.68 


0.90 


1 


SDSS J024933.42-083454.4 


46.38 


5230 


6601 


9.66 


0.25 


1.07 


4 


[HB89] 0329-385 


46.71 


7035 


8165 


10.11 


0.18 


0.86 


5 


[HB89] 0504+030 


46.32 


2046 


2402 


8.80 


1.52 


0.48 


6 


SDSS J100428.43+001825.6 


46.44 


3442 


3586 


9.33 


0.60 


1.53 


7 


TON 618 


47.31 


10527 


14905 


10.82 


0.14 


1.60 


8 


[HB89] 1246-057 


47.16 


5817 


6842 


10.22 


0.41 


1.82 


5 


[HB89] 1318-113 


46.89 


4150 


4898 


9.75 


0.61 


0.56 


5 


[HB89] 1346-036 


46.88 


5110 


5129 


9.94 


0.41 


0.78 


5 


SDSS J135445.66+002050.2 


46.49 


2627 


2208 


9.13 


1.08 


2.93 


7 


UM629 


46.56 


2621 


2521 


9.16 


1.15 


0.80 


7 


UM632 


46.54 


3614 


3828 


9.43 


0.60 


0.28 


6 


UM642 


46.29 


3925 


5478 


9.35 


0.40 


0.51 


7 


UM645 


46.31 


3966 


4972 


9.37 


0.40 


0.29 


7 


SBS 1425+606 


47.38 


3144 


6769 


9.82 


1.71 


0.47 


2 


SDSS J170102.18+612301.0 


46.34 


5760 


5944 


9.72 


0.20 


1.19 


7 


SDSS J173352.22+540030.5 


47.00 


3078 


5460 


9.57 


1.26 


1.79 


7 


[HB89] 2126-158 


47.25 


3078 


3354 


9.72 


1.58 


0.77 


2 


[HB89] 2132+014 


45.77 


2505 


2801 


8.65 


0.61 


0.80 


2 


2QZJ221814.4-300306 


46.54 


2986 


3070 


9.27 


0.88 


0.96 


1 


2QZ J222006.7-280324 


47.22 


5238 


5660 


10.16 


0.53 


1.67 


1 


[HB89] 2254+024 


46.45 


2597 


1954 


9.09 


1.07 


0.59 


6 


2QZJ231456.8-280102 


46.31 


3459 


2011 


9.26 


0.53 


1.21 


1 


2QZJ234510.3-293155 


46.32 


3908 


3453 


9.37 


0.42 


1.59 


1 



References. — 1. 2QZ archive (Croom et al. 2004). 2. Dietrich & Wilhelm-Erkens (2000). 3. Forster et al. (2001). 4. SDSS first data release (Abazajian et 
al. 2003). 5. Osmer & Smith (1977). 6. Baldwin, Wampler, & Gaskell (1989). 7. SDSS early data release (Stoughton et al. 2002). 8. Baldwin & Netzer (1978). 



line region (NLR), and a broad BLR component, that we will 
use to determine Mbh- Resolving each of those two physical 
components can be done if there is an obvious NLR spike in 
the H/3 profile, as is the case for several broad-line Seyfert 
1 galaxies (BLSls). In cases where a narrow spike is un- 
detected, the NLR contribution to the H/3 emission can in 
principal be estimated by taking FWHM([0 III] A5007) as its 
width, and its flux as a certain fraction of the [O III] A5007 
flux. This method is highly model-dependent, and in many 
of our sources we do not detect [O III] emission at all (see 
Paper II). The two Gaussians we fit to H/3 do not attempt to 
recover the physical NLR and BLR components, but merely 
to represent a fit to the entire profile. In luminous AGNs the 
NLR-to-BLR flux ratio, at least in the UV lines, is of the or- 
der of a few percent (e.g.. Wills et al. 1993), resulting in a 
negligible NLR contribution to the line flux. Since we also do 
not detect a significant NLR spike in any of our H/3 profiles, 
we assume that the observed line represents entirely the BLR 
contribution. 

We have also obtained another measure of line-width, based 
on the entire emission line flux, regardless of the profile. This 
is the inter-percentile velocity (IPV) width, which we have 
used for measuring the difference between wavelengths cor- 
responding to 25% and 75% of the integrated emission line 
flux across the entire line profile, starting with the continuum 
on the blue side and ending with the continuum on the red side 
of the line. Both the FWHM and IPV methods have their own 
merits, although it is still not clear which gives the closest 
estimate of the BLR velocity width (see Corbett et al. 2003 
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Fig. 4. — Monochromatic luminosities at 5100A vs. 1450A for the 
sample of 92 AGNs. Empty circles are broad-line AGNs (BLAGNs) with 
FWHM(H/3)>1500 km s"', and plus signs indicate 2 < z < 3.5 quasars. 
NLSls (with FWHM(H/3)<1500 km s"'; see SN02) are marked with (filled 
circles), and RLQs are marked with empty squares. The dashed line 
marks a 1 : 1 correspondence, and the solid line marks an empirical ratio of 
ALa(1450)/ALa(5100) = 1/(5100/1450) assuming f^ « v'" with a = 0.5. 
Most AGNs lie closer to the dashed line, indicating somewhat larger a. 



and references therein). In most cases, the IPV widths of 
H/3 are somewhat narrower than the corresponding FWHM 
(by ^15%; for a Gaussian profile the IPV width should be 
^43% narrower than the corresponding FWHM), although 
both methods give consistent line-widths within the uncertain- 
ties discussed above. 
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3.2. Luminosity, BH mass, and accretion rate 

We have derived monochromatic luminosities for our entire 
AGN sample at the rest-frame 5100A band using our new IR 
observations, and data from the literature, and compared them 
with the monochromatic luminosities at the rest-frame 1450A 
band obtained from the literature, and from the archives of the 
SDSS, and the 2dF quasar redshift survey (2QZ; Croom et al. 
2004). For ten sources we could not find 5100A fluxes, in- 
cluding six sources from our high-z sample (see Table^. In 
those cases we obtained the 5100A luminosity by extrapolat- 
ing from 1450A fluxes and assuming (as explained in §|3 a 
continuum function of the form: /^ <=<: i'~'^-^. Figure |4] shows 
this comparison, along with two lines, indicating a 1:1 cor- 
respondence, and a 1 : y^5100/1450 representing the expected 
ratio between the 1450A and the 5100A bands for a f^ oc i^~°' 
continuum with a = 0.5. Most AGNs lie in between the two 
lines indicating perhaps a somewhat redder {a > 0.5) contin- 
uum for our sample. We also note that most RLQs in our 
sample are slightly redder than the radio quiet objects, consis- 
tent with the steep continuum slopes in RLQs found by Netzer 
et al. (1995). Since the luminosities in both bands are highly 
correlated, and in order to avoid complications due to intrinsic 
reddening, we prefer to use the 5100A luminosities^ over the 
1450A luminosities in determining BH masses and accretion 
rates in our entire sample. 

Black hole masses and accretion rates were determined 
by the single-epoch mass determination method (Vestergaard 
2002), which is based on the Kaspi et al. (2000, hereafter 
KOO) /?BLR- ALa(5100) correlation, where ALa(5100) is the 
monochromatic luminosity at 5 100 A. In this method the BLR 
size is estimated from the measured ALa(5100) and the mass 
estimates follows from assuming Keplerian motion of the H/3 
emitting gas. The scaling used here is based on a re-analysis 

* For a few nearby Seyfert 1 galaxies in our sample, this choice might add 
contamination from the host galaxy. However, in most cases such a contribu- 
tion is much less than 50%, which is an acceptable uncertainty level, given 
the typical large expected continuum variations. 



of the KOO measurements (Peterson et al. 2004; Kaspi et al., in 
preparation) that gives ^blr °^ ^Lx (5 lOO)*''^ with a somewhat 
different size scaling factor. It also uses a somewhat different 
mass scaling factor that multiplies the FWHM(H/3) term in 
the mass determination equation (see KOO). The result. 



Mbh = 6.2 x lO*' 



ALa(5100) 
10^4 erg s-i 



0.6 r 



FWHM(H/3) 



n2 



103 km s- 



Mp 



(1) 
gives slightly larger masses compared with KOO and Netzer 

(2003) ^. We then convert the observed luminosity and the 
deduced masses to accretion rates, represented by LBoi/^Edd, 
by assuming LboI = 7 x ALa(5100). The factor 7 is in the 
middle of the range (5-9) representing observed and theoreti- 
cally deduced spectral energy distributions (SEDs; see Netzer 
2003). The result is 



ifiol/iEdd = 0.75 



ALa(5100) 



1044 



erg s 



0.4 



FWHM(H/3) 
10^ km s-i 



(2) 



Table |2] lists luminosity, Mbh, and iBoi/^Edd (hereafter 
i'/i'Edd) for our new 2 < z < 3.5 quasars. By combining the 
measurement uncertainti es o n FWHM(H/3), as given in Ta- 
bleland discussed in § 13.11 with the small (^ 25%) uncer- 
tainty on the luminosity, we obtain typical uncertainties on 
Mbh and on L/L-eaa, which are not larger than a factor of two, 
and therefore do not affect our main results as outlined below. 
A similar determination of Mbh, and L/Lndd was applied to 
LBQS 0256-0019, LBQS 0302-0000, and 61 sources from 
the SN02 sample. It is worth noting at this point that, within 
our new (albeit small) sample of luminous high-z quasars, 
there is no clear relation between Mbh and radio loudness (see 
Tables [n& |2}, as was previously suggested in several stud- 
ies, e.g., Franceschini, Vercellone, & Fabian (1998), and it is 
therefore consistent with the results of Woo & Urry (2002). 

' In fact, the most recent result of Kaspi et al., in prep., shows that the 
mass scaling factor used here is ^ 25% larger than theirs, but both factors are 
consistent within the uncertainties. 
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Fig. 6. — N V/C IV vs. accretion rate. Symbols are similar to those in Fig.|4] A strong metallicity-accretion rate correlation for all AGNs is apparent, and 
most NLSls are found in the same region of parameter space that is shared by the high-z quasars. Note that the location of Mrk 766 (the NLSl with the lowest 
accretion rate in this diagram) may be affected by strong intrinsic reddening (see text). 



3.3. BLR metallicity 

Following HF93, we use N V/C IV as the BLR metallicity 
indicator for our sample. This line ratio is easier to measure 
than most other emission lines in a typical spectrum, and was 
considered by Hamann et al. (2002) as a robust indicator since 
it minimizes non-abundance effects, such as temperature, and 
SED. About half of our quasars have rest-frame UV data 
available in the literature, and we complemented these data 
with measurements of N V and C IV for the rest of the sam- 
ple, which are mostly recently detected quasars. Rest-frame 
UV spectra of those sources were obtained from the SDSS and 
2QZ archives. The standard SDSS data allows reasonably ac- 
curate flux calibration and SED determination. The shape of 
the non-flux calibrated 2QZ spectra is more problematic and 
was obtained by dividing each spectrum by a quadratic, which 
provides a reasonable estimate to the spectral response func- 
tion of the 2dF (Folkes et al. 1999). Each SDSS and 2QZ 
spectrum was fitted with a power-law continuum, and two 
Gaussians for the C IV emission line. The C IV profile was 
then used as a template for all other rest-frame UV emission 
lines, and in particular N V. We list N V/C IVmeasurements 
for our 29 2 < z < 3.5 quasars in Table |2] The uncertainties 
on the measurement of this ratio vary between ~ 10% for the 
high-quality spectra and ~ 30% for the poorer quality ones 
(see also SN02). 

4. RESULTS 

4. 1 . A new metallicity-accretion rate relationship in AGNs 

In Figure|5lwe plot N V/C IV vs. (a) ALa(5100), and (b) 
Mbh for our sample of 92 AGN. One can see that the Z-L, 
and Z-Mbh relations show considerable scatter, mainly due 
to the NLSls which have low BH masses and luminosities. 
The first of those relations has been shown and discussed 
by SN02 for a slightly larger sample. The scatter is signif- 
icantly reduced once metallicity is plotted against accretion 
rate as shown in Fig. |6l In this case, all AGNs follow a 



Table 3. Spearman Rank Correlation Coefficients 



N v/C IV 


Number of 


'"S 


P 


vs. 


Sources 






(1) 


(2) 


(3) 


(4) 


L/L^i' 


92 


0.57 


< 10-5 


L/LEii' 


82'' 


0.54 


<io-5 


Mbh" 


92 


0.35 


3.1 X 10-^ 


Mbh" 


82'' 


0.53 


< 10-5 


ALa(5100) 


92 


0.53 


< 10-5 


ALa(5100) 


82'' 


0.71 


< 10-5 


L/Lem' 


82 


0.32 


1.6 X 10-3 


L/Lem'^ 


72'' 


0.41 


1.7 X 10-^ 



"Calculated from ALa(5100) and FWHM(H/3) 

''NLSls excluded 

"Calculated from ALa(1450) and FWHM(C IV) 



Strong metallicity-accretion rate (N V/C IV -L/Ledd; here- 
after Z-L/Ledd) correlation with NLSls being no exception. 

Table|3llists the Spearman rank correlation coefficients (rs) 
and the chance probabilities (P) for each of the three rela- 
tions. We also tested the correlations by removing the 10 
NLSls from the sample. As can be seen from Table |3j this 
increases rs considerably in the Z-L and Z-Mbh cases, but 
their removal from the Z-L/L^i^ diagram leaves rs almost 
unchanged. We note that formally, both the Z-L and the 
Z-Mbh correlations are significant (P <C 0.01). This rep- 
resents the fact that our sample contains only 10 NLSls, out 
of 92 sources, (representing well their fraction in the AGN 
population). However, this sub-group clearly lies away from 
the general trend in Figures [Sj^ & [Sj? emphasizing the fact 
that they do not share the same properties with other low 
luminosity/MsH AGNs. 

To summarize, our results indicate that the N V/C IV line 
ratio is correlated primarily with L/Ledd and not with lumi- 
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nosity or BH mass. It is the combination of luminosity and 
line-width (i.e., the accretion rate) that makes all the AGNs in 
our sample, including NLSls, follow the same trend. We cau- 
tion that our entire sample of 92 AGNs is neither complete 
nor fully representative of the AGN population as a whole, 
and that our main result may be subject to selection biases 
over the wide range of AGN properties, that were not fully 
explored in this work. 

4.2. FWHM(Hl3) vs. FWHM(C IV) as accretion rate 
indicators 

Reverberation BH masses were so far determined for 35 
AGNs using H/3 and the monochromatic luminosity at 5100A 
(KOO, Peterson et al. 2004), to set the empirical /?- L relation 
on which all Mbh estimates are made. For three AGNs, rever- 
beration masses were also determined using C IV (Peterson & 
Wandel 2000), and were consistent with the masses obtained 
from H/3. Vestergaard (2002) suggested that C IV can be used 
as a 'proxy' to H/3 and replace it as the reverberation Mbh 
indicator It is expected that the BLR components of both 
lines will satisfy FWHM(C IV)/FWHM(H/3) ~ ^/2, based on 
reverberation studies of both lines that show a smaller mean 
emission distance for C IV (Peterson & Wandel 2000). 

We have measurements of both H/3 and C IV in 82 AGNs 
of our sample, and the width of one line vs. the width of the 
other is plotted in Fig.0 One can clearly see that the widths 
of both lines do not follow a 1:1 correspondence nor a l:\/2 
dependence. In particular, NLSls with narrow H/3 lines show 
FWHM(C IV) which is not different than those of 'normal' 
BLSls, a fact that was first noted by Rodriguez-Pascual et 
al. (1997). Broad C IV in NLSls, and in other AGNs, has 
been attributed to high velocity gas outflowing from the center 
rendering a broad and blueshifted tail to this line (e.g., Leighly 
2000). 

To demonstrate the influence of FWHM(C IV) on NLSl 
Mbh estimates, we compared the Vestergaard (2002) single- 
epoch Mbh estimated from the UV spectrum with the re- 
verberation Mbh from KOO, directly measured from the H/3 
width, for four NLSls. We find that, on average, the Mbh 
UV estimates are larger by a factor of ~ 3 than the Mbh mea- 
sured by optical reverberation mapping. Hence, Mbh calcu- 
lated from C IV, at least for NLSls, may be systematically 
higher than Mbh estimates using H/3. An opposite behavior is 



also apparent in Fig.Q namely narrow C IV for broad H/3. 

To further illustrate the problematic use of C IV for individ- 
ual sources, we recalculated the BH mass and accretion rate, 
this time using the width of C IV, together with the Vester- 
gaard (2002) /?BLR-i relation. Figure |8] plots N V/C IV vs. 
accretion rate, and Table |3] lists the corresponding Spearman 
rank correlation coefficients. By inspection of Fig.|8land Ta- 
ble |3]it is apparent that the correlation between the accretion 
rate calculated from the C IV width, and N V/C IV, is much 
weaker compared with the accretion rate calculated from H/3, 
and that the Z-L/Leaa now has much more scatter. This re- 
sult may also explain the lack of a Z-L/Lem correlation in 
Warner, Hamann, & Dietiich (2004). 

5. DISCUSSION 

The main result presented in this paper is a strong correla- 
tion between metallicity (a la HF93) and accretion rate in a 
sample of 92 AGNs spanning a broad range of AGN luminos- 
ity (~ lO'*^- lO'*^ erg s~'), line-width, and radio loudness. We 
also confirm the earlier SN02 finding that AGN samples that 
include NLSls do not show a clear Z-L correlation. The new 
relationship, although supported by a strong correlation, is not 
free of scatter, and the effects of other physical properties are 
still unknown. In this section we discuss potential sources for 
intrinsic scatter in the Z-L/Ledd diagram as well as poten- 
tial caveats in the derivation of BLR metallicity and accretion 
rate. We also comment on the potential implications of our 
new result on metal enrichment and BH growth scenarios in 
galactic nuclei, and on the starburst-AGN connection. 

5.1. How accurate is the BH mass determination at high 
luminosity? 

Our estimate of BH mass, and hence accretion rate, rely on 
the KOO relationship applied to a sample of 34 AGNs with 
BH masses measured via reverberation mapping and modi- 
fied, slightly, in Kaspi et al., in preparation. That relation- 
ship is limited by the small number of sources, the luminosity 
range (^ 10^^- 10^^), the assumed constant SED, and perhaps 
other AGN properties. Mass estimates for sources beyond 
this luminosity range are therefore uncertain. Such extrap- 
olations lead to BH masses exceeding 10"^M3 in this work 
and in Netzer (2003). Hence, if we have somewhat overesti- 
mated Mbh in our luminous AGN sample, then the accretion 
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rate must have been underestimated. This will move lumi- 
nous quasars to the right hand side of Fig.|6land will produce 
an even stronger Z-L/Lgcid correlation. However, a very dif- 
ferent L-Mbh relationship at high luminosity may introduce 
more scatter into some of the correlations presented here. 

5.2. Is N V/C IV an adequate metallicity indicator for all 
AGNs? 

Using a grid of photoionization calculations, performed 
with the locally optimally emitting cloud model (LOC; e.g., 
Baldwin et al. 1995), Hamann et al. (2002) found that 
N V/C IV can trace the BLR gas metallicity, and it only weakly 
depends on the chosen SED and density. Our new results 
show that N V/C IV depends primarily on L/Lgcid, but it is not 
yet clear whether high N V/C IV is indeed indicative of high 
metallicities. In particular, it is crucial to compare the obser- 
vations of this line ratio to those of the model-independent 
weak-lines proposed by Shields (1976). Our investigation of 
this issue is beyond the scope of this work (see also Dietrich 
et al. 2003). 

Part of the scatter in the Z - L/L^dd diagram may be con- 
tributed by the yet unknown effect of BLR density on the way 
BLR metallicity is derived. The problem has been mentioned 
by HF93, but the observational verification is still missing. It 
is reasonable to assume that high metal contents are associ- 
ated with high densities in the nuclear region which are indi- 
cators of recent star formation. In particular, NLSls which 
have high ratios of N V/C IV, are also known to have high 
BLR densities, derived from the Si III] A1892/C III] A 1909 
density-sensitive line ratio (e.g.. Wills et al. 1999). We have 
measured this line ratio in our sample (available for about half 
of our sources), and found no clear trend with either metallic- 
ity or accretion rate. The two emission lines used for this ratio, 
which are fairly close to each other, are hard to measure and 
to deblend. Therefore the potential effect of the BLR density 
remains unclear at this stage. 

The continuum SED may also affect the way N V/C IV de- 
pends on BLR metallicity. It is not yet clear whether the SED 
depends on luminosity, or perhaps the accretion rate, but if 
the SED in high accretion rate sources is significantly differ- 
ent from the ones used by Hamann et al. (2002) to calculate 
N V/C IV, then it may cause some deviation from the canon- 
ical Z-L dependence. In particular, NLSls appear to have 
SEDs which are somewhat different from 'standard' AGN 
SEDs (e.g., Romano et al. 2004), and related perhaps to the 
putative (hotter) slim accretion disk (e.g., Wang & Netzer 
2003). As explained in § 13.21 a different SED may also af- 
fect our derived L/Ledd- 

A viable model of AGN evolution is that they start off with 
high accretion rates, when fuel supply is abundant, and grad- 
ually move towards lower accretion rates as the fuel supply 
dwindles. If the metal content of the gas does not change as 
a consequence of the fuel shortage, one can expect an AGN 
to move horizontally, from right to left, across the Z-L/L^dd 
diagram (Fig.|5), as a function of time. In other words, similar 
metallicities may be observed for a range of accretion rates. 
The dependence on age may be further complicated by the 
AGN duty cycle (assumed to be ^ 10^-10^ yr), as we do not 
know the number of cycles each AGN in our sample has gone 
through. If metal enrichment is strongly related to AGN fuel- 
ing, then a number of active episodes may contribute to metal 
enrichment provided the metal rich gas is not being depleted 
or blown out of the nucleus. This can lead to a Z-L/Ledd 
dependence which is more compUcated than assumed here. 



5.3. Intrinsic reddening 

Intrinsic reddening may affect the location of an AGN on 
the Z-L/Ledd diagram, by moving it diagonally towards the 
lower left corner in Fig.|6l due to weaker N V with respect to 
C IV and lower luminosity. Correcting for intrinsic reddening 
is problematic, and requires a broader spectral range to com- 
pare suitable line ratios (e.g., Crenshaw et al. 2002). We did 
not attempt to correct the data for intrinsic reddening but note 
that in the case of NLSls, such a correction may be necessary, 
since these objects are suspected to possess large amounts of 
dust and gas (e.g., Crenshaw, Kraemer, & Gabel 2003). We 
note that at least one NLS 1 in our sample (Mrk 766) may be 
strongly affected by intrinsic reddening, which might explain 
its location in Fig.|6l away from the rest of the NLSl popula- 
tion. 

5.4. Metal enrichment and BH growth in galactic centers 

If N V/C IV is an adequate BLR metallicity indicator, then 
our Z-L/Ledd correlation may have important implications 
for the understanding of BH growth in galactic nuclei as well 
as for the starburst-AGN connection. We find that luminous 
quasars at high-redshift, as well as NLSls in the nearby Uni- 
verse, may be experiencing (or have recently experienced) a 
vigorous star-forming episode, and their BHs are in the stage 
of rapid growth. In this sense, the NLSls in our sample may 
be in the early stages of their (current) activity, i.e., relatively 
'young' systems, as suggested by Mathur (2000). All the 
high-z quasars in our sample also exhibit high accretion rates 
and high metallicities and may be in a similar evolutionary 
stage. These sources are thus the high luminosity analogs of 
NLS 1 s and should, perhaps, be referred to as narrow-line type 
1 quasars (NLQls). However, our sample does not contain 
faint quasars at large distances and those must be added to the 
sample to test this idea. In particular, our new Z-L/Lgdd rela- 
tionship may be affected by a potential selection effect against 
low accretion rate AGNs at high-z, if such sources are found 
to be metal-rich. Our new result is consistent with a scenario 
in which there is an intimate relation between AGN fueling, 
manifested by accretion rate, and starburst, responsible for the 
metal-enrichment. An obvious avenue to proceed is to search 
for other evidence of recent star formation in our high-z sam- 
ple using IR, mm, and sub-mm observations. 

6. CONCLUSIONS 

We present new near-lR spectroscopic measurements of the 
H/3 region for a sample of 29 luminous 2 < z < 3.5 quasars. 
This sample is augmented with measurements of other AGNs 
(92 AGNs in total) for which spectroscopic data of the rest- 
frame UV band as well as the H/3 region are available. We 
show that BLR metallicity, indicated by N V/C IV, is primar- 
ily correlated with the accretion rate, indicated by luminosity 
and FWHM(H/3). We also show that using C IV as a proxy 
to H/3 may be problematic and can lead to spurious BH mass 
estimates. The potential implications of our new result to the 
starburst-AGN connection and to AGN fueling are briefly dis- 
cussed. 
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